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ABSTRACT 


A shock wave passing through the solar wind can apparently 
reflect energetic protons to build up an enhanced Intensity in 
front of the shock wave. Values of the difijsion coefficient in 
interplanetary space and the reflection efficiency of the par- 
ticles can respectively be calculated front the distance the 
enhanced region extends in front of the shock wave and the in- 
crease in intensity fron the undisturbed intensity. Values of 
2 

^ 10^7 £S_ and 90X are found by fitting the date of Ogilvie 
sec 

and Arens (1970) and describe properties of I MeV protons with 
pitch angle fn9(P scattering throu^ 20^. 


INTRODLCTION 


Singer (1970), Ogilvie and Arens (1970), and Armstrong, Krimigis and 
Bchannou (1970) have seen enhanced energetic proton Intensities for 15- 
20 minutes before the passage of a shock wave passing through the solar 
wind. Axford and Reid (1963) proposed a Fermi acceleration mechanism 
using the interplanetary and earth's bow shock waves to produce such 
enhanced intensities. 

Fermi acceleration between the interplanetary shock wave and scattering 
centers in the solar wind ahead of the shock wave (Van Allen and Ness, 

1967; Fisk, 1970) could also produce such enhancements. Here wo use ran- 
Qom walk theory in the presence of the reflecting shock wave and scattering 
centers to calculate curves of enhanced intensity vs. distance in front 
of the shock wave. The scattering centers are magnetic field irregularities. 

First we will ignore the directionality of the magnetic field and 
look at the enhanced region in the same way as one would an aerosol cloud 
in front of a moving screen which lets air molecules through but is 
impermeable to aerosol particles. We consider the case where the aerosol 
particles have speeds much larger than the screen speed. The particles 
bounce off the screen and random walk due to collisions with air molecules. 
Most particles will not move far from the screen because the net mean 
velocity (net distance traveled /time) of random walking (diffusing) 
particles decreases with time, so the screen will eventually catch up 
with each particle and reflect it again. So we will have a cloud in 


front of the screen which will not diminish with time. 


Next, in this paper, ve will consider reflection efficiencies of 1 
MeV protons off shock waves and show how increased or decreased intensity 
regions can exist in front of shock waves. 

Thirdly, we will talk briefly about the effects of the directionality 
of the magnetic field. Last we compare the calculations with experimental 
data and find a value for the diffusion coefficient of 1 MeV protons at 
1 A.U. 

calculations 

\ 

The enhanced Intantity region in front of the shock wave should be 
time invariant since the 15*20 minute interval needed for the passage of 
the enhanced intensity region is much smaller *han the time taken by the 
shock wave to travel from the vicinity of the sun to the satellite, near 
the earth. Buildup and loss rates are hypothesised to be equal because 
the peak intensity is only several or several tens of times hl^er than 
the background intensity that existed before the buildup. If the build- 
up race were higher than the loss rate, we mi^t expect a peak to back- 
ground ratio of at least several times 

^ z i so 

since energisation by elastic scattering off the mov: ig shock wave would 
considerably enlarge ti^is nun^er. Here the approximation is made that 
Che shock wave is planar, ell particles are piled up, and 

d^ ■ distaneo shock wave has traveled throu^ solar wind, 
d 2 ■ distance enhanced region extends in front of shock wave. 
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\ 

V = shock wave speed through solar wind ~ 200 

v„ ® solar wind spef;d ~ 400 ^ , 

=» - sec 

S 

L * distance to vicinity of sun ^ 10 km, 

tj^ = time taken for enhanced region to pass by satellite ~ 1000 sec. 

The intensity in front of the shock wave, Ig. should depend only on the 
distance from the shock wave, x, making ls(x) a function of one dimension. 

Diffusion Coefficient 

The steady state solution in front of the shock wave should roughly 
approximate the distribution obtained by a random walk process from a 
stationary point after a time T such that the meau random walk distance 
equals the distance traveled by the shock wave. The random walk distri- 
bution and mean distance are 



0 ' a, 
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is the diffusion coefficient. is the mean free path and V is the speed 
of a particle. We set 



Figure 1 shows how good an approximation this curve is compared to the 
more exact distribution, \^(x,0), c Iculated below: 

Since the intensity in front of the shock wave, Ig(x), is unchanging 
with time , we want a solution such that 
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where the shock wave moves a distance 6 in time t as shown in Figure 2. 
We let all particles random walk with perfect reflection off the shock 
wave to conserve the total number of particles since we first assume the 
pre-shock undisturbed interplanetary intensity, Ij,» is zero. Later, 
when we calculate the reflection efficiency, we will make the loss by 
absorption into the shock wave equal to the gain by picking up new 
particles from the undisturbed intensity and by energization from re- 
flection off the moving shock wave in order to conserve the total number 
of particles. We now let the shock wave jump a small distance 6 at 
time t“0 and sweep up all the particles in the space traversed. These 
particles then random walk from x=6 during the time from t*0 to t=T. 
Meanwhile particles at x > 6 random walk in the presence of a reflecting 
barrier at x=6. Using random, walk theory. 
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The bracketed factor In the first Integral gives the distribution for 
random walk with a reflecting barrier (Chandrasekhar, 1943). The second 
integral is the sou-ce of parti. les that random walk from x«6. 

The approximation of the shock wave jumping and putting the source 
at x>6 should give the same result as a smoothly progressing shock wave 
because the random walk motion of particles is much faster than the shock 
wave speed at first; i.e., 


'-N — 

X 

7't 
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where 


i i r,„,i 




c y 


and is the time in which the average number of collisions each 

particle makes with the shock wave is 1 as shown in Figure 3. If the 
solution for y(x,0) satisfies 


0 i ( 
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this approximation will be valid. 
We approximate 


^(vVS, r) 
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and Integrate by parts: 



if we neglect terms of order 5^. 


Iterative solutions to V(y 0) might be expected to converge if an 
initial distribution that goes to ?.ero for large x slowly evolves to one 
final distr ib’.cion as the shock wave progressrr. This fact is borne out 
by solutions obtained by iteration on a computer and shown in Figure 1. 

j also seen to be finite, justifying the approximations 

made near the origin. 

Reflection Efficiency 

Since we hypothesize a steady state solution, we require that the 
total number of narticles in the enhanced region be a constant. If the 
intensity in the enhanced region is 

IsiE) “ 

where E is energy and |j, and y are constants and each particle has its 
energy multiplied k times upon reflecting off the shock wave, the intens 
would be increased by k^"^ after all particles have been reflected one 
time. If € is the efficiency for reflection, we need 

e = 1 

for no net increase in intensity since we want a steady state solution. 
In general e will not satisfy this equation; define to such th^t 

k>"^ Co « 1. 

We will first work with in order co gain understanding of the problem 
and then find a solution for e. 



- 9 - 

The equation for is not entirely correct if there is a background, 
before the enhanced region arrives. The shock wave must let as 
many protons pass througf. it as it both sweets up and creates by energi- 
zation if the enhanced region is to remain unchanging. If I 3 .. (E) is the 
intensity per unit energy interval that is moving toward the shock wave 
just in front of it as seen in the solar wind rest frame and t <'<' 

+ v) T is the number per unit energy interval that hit the 
shock wave in time t. Vj^ and Vg are the average velocities of the par- 
ticle” comprising I 5 . perpendicular to and parallel to the shoes wave 
front. 

Approximation, Is_(E)(V^^ + v) TC^k is the intensity per unit 
energy interval recc'ding from Che shock wave after striking it in time 
If the uumber of particles is to remain constant. 

If 

and Cq is large, then 

Js-U) - i ^sCe) 

Hudson (I965) and Sonnerup (1969) have devised models that describe 
particles re..lecCing off shock waves. Here we will use a simple approach 


and 
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and find k as follows: For elastic reflection and isotropic fluxes 
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if 







and where superscript p denotes post ccHision.- 

As an example, the velocities of a 1 MeV proton and a typical shock 
wave are 


V U ti*v- 

" I.-5S 


"LT r: .^OO 


yielding 


If 






J- ^ 1 


\o 





Such a large value for Cg can seem plausible for both thick and thin 
shock waves. For a thick wave draw an analogy with grazing incidence 
energetic electrons scattering off a mefal surface, where reflection 
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efficiencies ~ 997, can often occur. An electron will scatter through a 
small angle. If it scatters out of the block it is free; if ir scatters 
in it still has a large probability of subsequently scattering out before 
losing all its energy. If the magnetic field inside the shock wave in 
our problem is much more turbulent than in the solar wind ahead of it, 
the protons will strongly scatter inside the wave until they become "free” 
by exiting from the wave. These protons, of course, do not gain or lose 
much energy and hence there is practically no limit on the time they can 
stay in the wave. If the thickness of the shock wave is larger than the 
mean free path, 1, of a proton, the protons will predominately exit on 
the side they entered; i.e., the front of the shock wave. 

If the shock wuve is much thinner t'nan the gyroradius of a particle, 
the particles will hardly notice the shock wave and have a very low re- 
flection coefficient. However, almost as many particles will overtake 
the shock wave and pass through the wave from the beck to the front side, 
making the effective reflection coefficient < 1 if the intensity on the 
back side is equal to the undisturbed intensity on the front side. But 
now we cannot calculate an effective energization factor. More precise 
calculations considering magnetic fields will be considered later. We 
will proceed using the thick shock wave model. 

Returning to the actual reflection coefficient, e must fall within 
certain limits in order to give an unchanging region in front of the 
shock wave. If 

^ ^ C aVa \ I 
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there will be a lower intensity in front of the shock w-ve than in the 
background; i.e., the sh 'k wave acts as a sink for particles. The 
decrease seen by Van Allen and Ness (1967) may have been due to an 
absorbing shock wave. If 

5 > k 

a steady state solution cannot exist. Particles are created by energization 
faster than they are lost into the shock wave. The enhanced region will 
constantly increase as the shock wax'e progresses. So e must satisfy 



for a steady state solution with 

Is(E) > Ib(E). 

e most certainly almost never falls within these close limits and 
hence steady state enhanced regions in front of shock waves as here des- 
cribed should seldom exist. Several possibilities might explain the 
data: 

1) The enhanced intensity is not unvarying but is Increasing and 
the reason it is not extremely large is that e or D changes 
with distance from the sun so that the enhanced intensity has 
been increasing for only a short while. 

2) Particles leak around the edge of the shock wave. 

3) The shock wave surface is irregular with regions that act as 
sinks (low e) to moderate the growth of enhanced intensities 
in other regions (high ?). 
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U) the <\xford-Reid forrai acceleration mechanism is the only mechanism 
that causes enhanced intensities in front of shock waves. 

5) There is a finite source of low energy particles caused by de- 
creasing e with decreasing energy. 

Let us consider the fifth possibility. We suppose that below some 
energy, E©, the shock wave does not reflect and energize particles. If 
interplanetary space were devoid of particles with energy > E© except for 
one narrow band of energies of width 

ae'- (Uf)-l) -=-+ fc -- f ' 

where the intensity (per unit energy per unit distance) is IjjCE'), the 
shock wave would see a steady source of particles as it swept through 
space. These particles would be energized by the shock wave to produce 
a distribution in front of the wave. AE ' was chosen so that after the 
first reflection the new intensity would fall adjacent to the original 
intensity as shown in Figure 4. 

The fraction of particles in the enhanced region at any energy 
striking the shock wave in unit ti-ie is 

TO i T, Lt,o) V / f r, Cr, J X (i) 

Cell ' c, 

where we also consider the shape of the enhanced region. For an unchanging 
enhanced region in front of the shock wave , each of the blocks in Figure 
4 must receive particles from the next lower block at the same rate as 
it loses to the next higher block: 

Cc^ o) IS r ~ < . 
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Using the analytic approxii.iate expression P(x,l) for the enhanced intensity 
distribution Ig(E,x) we get from (1) 



In the time the average number of collisions each particle makes 

with the shock wave is 1. In this time we energize all particles with 
energy E in the enhanced region once on the average. In the time 
we admit I{,(E') v Tj.qJ^j^(E') particles from the undisturbed region into 
the enhanced region. So all particles move to the right one block in 
Figure 4 and we obtain an approximate intensity at the shock wave by 
using P(x,T) 


U) 


T, (r, .y-w - >£ 


The factor containing k(E) accounts foe overlapping blocks after the 
first reflection in Figure 4. 



is the efficiency for particles starting at energy E* to be energized 
to E without being lost by absorption into the shock wave, where N is 
the number of reflections and e is assumed to be independent of energy. 
We can approximate R(E,E') for e 1 and k es 1: 
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Now, actually the source of '':i- .Icles is not at one narrow band of 
energies but ir- the undisturbed interplanetary intensity, 

■ 31, 

assuming a power law spectrum. We obtain 
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So we obtain 



How this result depends on the various variables is shown in Figure 5. 

MAGNETIC FIELD AND PITCH ANGLE 

We have neglected the magnetic field in the solar wind in all the 
above discussions. Actually the reflection efficiency depends on the 
angles the magnetic field lines in front of and behind the shock wave 
make with the shock wave and on the pitch angle and the phase of the 
particle when it hits the shock wave (Hudson, 1965). In the following 
analysis we will make some very rough approximations to find the general 
shape of the pitch angle distribution in the enhanced region. 
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If there is no correlation between the shock wave normal and the 
magnetic field direction, the average angle between the two is 1 radian. 
Looking at the computer calculated graphs of Hudson we see that for this 
angle, particles with pitch angle 70° <” v 110° are all reflected while 
those with smaller or greater pitch angle are less efficiently reflected. 
In an isotropic flux, particles in the above mentioned pitch angle range 
comprise about 1/3 of the total number of particles. If the intensity 
of these particles is enhanced at the shock by a factor of 20 over the 
undisturbed isotropic intensity and particles with other pitch angles are 
enhanced by a factor of 5, the enhanced intensity averaged over pitch 
angle would be increased a factor of 10 over background. The intensity 
of particles with pitch angle s« 90° would be several times the intensity 
of particles with velocity more aligned with the field lines. 

The above numbers have not been rigorously calculated; pitch angle 
data is not yet available making comparisons between theory and experi- 
ment impossible. When pitch angle measurements become available, a more 
exact analysis might be useful. However, the calculations performed in 
this paper along with Hudson's analysis (1965) predict the general fea- 
tures of the pitch angle distribution described here. 

COMPARISON WITH DATA 

Is(E,x) is very sensitive to a, e, and E/Eq which perhaps accounts 

for the wide range of peak heights we see in the data mentioned in the 

introduction. In the solar wind rest frame a 200 eV proton has roughly 

the same velocity as a 200 ~ shock wave, which would thus place an 

sec 
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approximate lower limit on Eq. The approximation for R(E, E') is very 
poor for large E/Eq, however it can be seen from Figure 5 that e must be 
in the neighborhood of 0.9 for 


Of = 3 



which are reasonable values for 1 MeV protons (Ogilvie and Arens, 1970). 
Using the lower x-scale in Figure 1, a plot of 'y(x,0) with 

1A 2 

D = 6.0 X 10^'^ SID- 
sec 

is shown superposed on data from the November 29 shock wave. The x-scale 

was obtained by knowing that the data points were taken at 82 second 

intervals, the shock wave speed v was 52 ^22- relative to the solar wind, 

sec 

the solar wind speed v„„ was 386 ISID., , and the angle between v and v-u was 

aw ggg Bw 

37° (Ogilvie and Burlaga, 1969. N.B. Correctio.iS have been made in this 
paper) . 

where D|| and are the parallel and perpendicular diffusion coefficients 
and Q is the angle between the magnetic field and the normal to the shock 
wave. ^ was for the November 29 event and if D|| is much larger than 





If th6 enhanced intensity is largely made up of particles with pitch 

angle 

80° < j < 100°, 

the mean free path along the field line is approximately 
/■) ^ ^ A ^ f 

-^(1 ' " 5 - 5*00 U ^> 

The gyroradius of a 1 MeV proton in a 5 ganina field is 3 x 10^ km. The 
result for the mean free path implies that a proton will scatter through 
an angle ets 20° once every gyrore volution. 

The mean free path and the diffusion coefficient found here are 
much smaller than the values found by considering the time needed for 
particles to propagate from a solar flare to the vicinity of the earth 
(Bryant et al, 1962; Hofmann et al, 1963; Krimigis, 1965; McCracken et al, 
1967). Our mean free path denotes a distance traveled along a field line 
by a particle near the earth with pitch angle 90° before it scatters 
through an angle 20°. The other determinations o. the mean fret path 
are averages of the mean free path between the si>n and the earth rather 
than the mean free path near the earth. Also, in those paper particles 
with velocities more aligned along the field lines are i,. eluded and the 
scattering angle is much larger, being 180°. 

Similar calculations for other shock waves are difficult to make 
because the events are complicated or, more often, some necessary mea- 


surements are lacking. 
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FIGirRE CAPTIONS 


1. Approximate (P(x,T)) and exact (Ig(x)) intensity distributions in 
front of the shock wave. Data from three shock waves (Ogilvie 

and Arens, 1970) are also plotted with the x-scale spac.ng adjusted 
to make reasonable agreement with the theoretical curves. Data 
points have the undisturbed intensity, I^, subtracted and are 
normalized to the origin. Data points were taken at 82 second 
intervals, which is the sampling rate of the satellite experiment. 

2. Unchanging intensity distribution, t(^)> ^ shock wave moves. 

3. Paths of shock wave and typical particle. 

A. Intensity distribution of a batch of particles after successive 
reflections off the shock wave. 

5. Intensity at the shock wave is very sensitive to a, e and E/E^. 
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FIGURE 5 



